Finite density QCD with heavy quarks 
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In the large fermion mass limit of QCD at finite density the structure of the partition function greatly simplifies 
and can be studied analytically. We show that, contrary to general wisdom, the phase of the Dirac determinant 
is relevant only at finite temperature and can be neglected for zero temperature fields. 
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Introduction 

Finite density QCD is still a poorly understood 
field from an analytical point of view. We miss 
any quantitative data about the passage from the 
hadronic to the quark-gluon plasma phase at high 
densities and even more qualitative information 
like the order of the phase transition (if any) is 
not available. 

This situation mainly reflects the difficulties 
to perform numerical simulations of QCD at fi- 
nite chemical potential /i. Almost all the stan- 
dard algorithms for zero density QCD include the 
fermionic determinant in the generation of con- 
figurations and fail for complex valued determi- 
nants. Attempts to create ad hoc algorithms have 
been made but they are generally very computer 
consuming, they rely on some a priori unjusti- 
fied approximations to get real valued observables 
and the results obtained up to now are plagued 
from unphysical effects like the early onset of the 
baryon density Q , M . 

The main difficulty is to evaluate correctly the 
contribution of the phase of the action. In fact the 
partition function can be written as the product 
of two terms: 



Z||(cOS 



(1) 



where the first is evaluated using the modulus of 
the Dirac determinant and the second is the mean 



value of the cosine of the phase calculated weight- 
ing the configurations with the usual pure gauge 
part of the action times the modulus of the de- 
terminant. In the thermodynamic limit (cos 0) | 
gives a finite contribution to the free energy den- 
sity if it is proportional to e~ v . If this is the case 
there is no chance to determine numerically its 
value: for any reasonable statistic, is not possible 
to evaluate it from a set of measurements of O(l). 

The simple idea to consider only the modulus 
of the action has received very little attention, 
especially after the studies on the random ma- 
trix model where the phase of the Dirac deter- 
minant is fundamental for obtaining the correct 
result ||. In QCD the one dimensional case can 
be addressed analytically and we can see that the 
average value of the cosine of the phase goes to 
1 as we approach the thermodynamic limit [Q; 
the physically relevant four dimensional theory 
can only be studied in the large bare mass limit 
where several simplifications in the fermionic de- 
terminant structure occur. 



The strong coupling limit 

We can first write the fermionic matrix A sep- 
arating the contribution of forward (backward) 
temporal links G (G T ) and spatial links V: 
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2A = 2ml + e^G + e^G* + V 



(2) 



When to 3> 1 only the mass term and the forward 
temporal links can contribute to the fermionic de- 
terminant. Since quarks do not propagate det A 
can be written as a function of traces of Polyakov 
loops (Li) only and we get, except for an irrele- 
vant multiplicative constant, 



det A = e 



3V S L 



c 2 TrLj 



+ cTrL*i +1)(3) 



where V s is the spatial volume, L t the lattice tem- 
poral extent and c = (|?) Lt §. 

In the strong coupling limit we can solve ex- 
actly the theory since the Polyakov lines are in- 
dependent and the partition function factorizes 
trivially: 

Z(fi = 0) = e 3V ° L ^ (c 3 + l) Vs (4) 

It follows immediately that, in the thermody- 
namic limit, the system undergoes a first order 
transition at /i c = log 2m where the baryon den- 
sity jumps from zero to the saturation value. 
Whatever the spatial volume is the singular be- 
haviour in the free energy density is related to L t 
only. 

Regarding equation ([!]) we see that an upper 
bound for Zu is given by the partition function at 
[3 = oo and therefore 
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As the volume goes to infinity this ratio goes to 
1 except when [i — [i c where it goes to zero expo- 
nentially with the spatial volume; we have a cor- 
rect description of our zero temperature system 
replacing the fermion determinant by its absolute 
value. 

The situation changes completely if we con- 
sider the effect of non zero temperature: we 
must keep L t finite and consider the limit of in- 
finite spatial volume. From (||) we see that the 
transition becomes a smooth crossover and, from 
the definition (|l|), we can also realize easily that 
(cos 0} || ~ e~~ Va . Since the system volume is pro- 
portional to V s this implies that the phase can 
not be neglected. 

The availability of the analytical solution ([|) 
can also provide a useful tool to check the re- 
sults of numerical simulations. Similarly to the 
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Figure 1. Logarithm of the modulus of the 
GCPF coefficients times the sign of the coeffi- 
cients as a function of the coefficient order: ana- 
lytic (crosses) and numerical (diamonds) results 
for V = 4 3 x 4. 



Grand Canonical Partition Function formalism 
we have written the partition function Z as a 
polynomial in the variable c and studied the con- 
vergence properties of the coefficients. We gen- 
erated several thousands of random gauge config- 
urations and calculated the averaged expansion 
coefficients. This has been done for lattices 4 4 
(fig. 1), 4 3 x 20 (fig. 2) and 10 3 x 4 (fig. 3) 
and results are plotted, as a function of the coef- 
ficient order, superimposed to the ones obtained 
from (Q). From the plots we can realize that i) 
numerical results are not positive definite as they 
should be but for large L t this phenomena disap- 
pears; ii) even considering the modulus of the av- 
eraged coefficients we do not get the correct order 
of magnitude in the central region; Hi) increasing 
L t or V s the situation makes worse. 

This seems to indicate that, in this model, 
our configuration ensemble is not relevant for the 
physics at intermediate c (i.e. c ~ 1) and a correct 
determination of all the coefficients from a ran- 
dom ensemble is not possible. Increasing the tem- 
poral lattice extent the contribution of the phase 
becomes smaller and, with reasonable statistic, 
we get positive averaged coefficients. Neverthe- 




Figure 2. As in fig. 1 for V = 4 3 x 20. 



Figure 3. As in fig. 1: analytic (uppermost 
continuous line) and numerical results for V = 
10 3 x 4. 



less their convergence to the correct values do not 
improve suggesting that a non accurate determi- 
nation of (cos 4>) || may not be the only source of 
artifacts in simulations (see also |Q). 

Finite coupling 

At non zero (3 the pure gauge part of the action 
couples the Polyakov lines among themselves and 
any trivial factorization in the partition function 
is not possible any more. The partition function 
can be written as the product of the pure gauge 
partition function times the (3 — partition func- 
tion times a factor R: 
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(6) 



From the previous section we conclude that, in 
the zero temperature case, the only possible con- 
tribution of the phase comes from @ . For large 
temporal extent c approaches or oo and in both 
cases the absolute value of all the factors in the 
numerator of (^J) is upper and lower bounded by 
a finite number. This implies that R behaves at 
most exponentially with the spatial lattice vol- 
ume V s except for fi — log 2m. When this hap- 
pens the phase is clearly not relevant and the 
free energy density reduces to the pure gauge 
contribution plus the infinite coupling term: we 



get again a first order saturation transition at 
/i c = log 2m for all the values of (3. 

At finite temperature the scenario is completely 
different. The non analyticity in the Z((3 = 0) 
disappears and expression ^ can also give a fi- 
nite contribution in the free energy density. Even 
if an exact evaluation of R is not possible this 
points toward a smooth crossover from the zero 
density to the saturation regime as soon as L t is 
considered finite in the thermodynamic limit. 
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